Non-Newtonian viscosity of E-coli suspensions 
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The viscosity of an active suspension of E-Coli bacteria is determined experimentally in the dilute 
and semi dilute regime using a Y shaped micro-fluidic channel. From the position of the interface 
between the pure suspending fluid and the suspension, we identify rheo-thickening and rheo-thinning 
regimes as well as situations at low shear rate where the viscosity of the bacteria suspension can be 
lower than the viscosity of the suspending fluid. In addition, bacteria concentration and velocity 
profiles in the bulk are directly measured in the micro-channel. 

PACS numbers: 47.63.-b, 47.57.Qk, 47.57.E- 



The fluid mechanics of microscopic swimmers in sus- 
pension have been widely studied in recent years. Bac- 
teria [3, 13], algae @, H[ or artificial swimmers Q dis- 
persed in a fluid display properties that differ strongly 
from those of passive suspensions [f|. The physical re- 
lationships governing momentum and energy transfer as 
well as constitutive equations vary drastically for these 
suspensions 0, § ■ Unique physical phenomena caused by 
the activity of swimmers were recently identified such as 
enhanced Brownian diffusivity [l], Isl— Tioj] uncommon vis- 
cosity [i , 12 , 13| , active transport and mixing fill or the 
extraction of work from isothermal fluctuations |l3l. [Tij . 
The presence of living and cooperative species may also 
induce collective motion and organization at the meso- 

impacting the consti- 



scopic or macroscopic level [13 Ul 
tutive relationships in the semi-diluted or dense regimes. 
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FIG. 1: Experimental set-up. (a) Time averaged image of 
the micro-channel (W — 600 um) for Q — 10 nL/sec and vol- 
ume fraction <f> = 0.35%. The red and blue frames indicate 
the measurement areas, (b) Concentration profile C(y) nor- 
malized by the maximal concentration Cm (black line) and 
error function fit used to determine the interface position (red 
line). 



The E.Coli bacterium possesses a quite sophisticated 
propulsion apparatus consisting of a collection of flagella 
(7-10 fim length) org anized in a bundle and rotating 
counter-clockwise pol ]. In a fluid at rest, the wild-type 
strain used here has the ability to change direction by 



unwinding some flagella and moving them in order to 
alter its swimming direction (a tumble) approximately 
once every second plj . In spite of the inherent com- 
plexity of the propulsion features, low Reynolds number 
hydrodynamics impose a long range flow field which can 
be modeled as an effective force dipole. Due to the thrust 
coming from the rear, E.coli are described as "pushers", 
hence defining a sign for the force dipole which has a 
crucial importance on the rheology of active suspensions 
31 . For a dilute suspension of force dipoles, Haines et al 
22j and Saintillan [24| derived an explicit relation relat- 
ing viscosity and shear rate. They obtained an effective 
viscosity similar in form to the classical Einstein relation 
for dilute suspensions : rj = 770 (1 + K<f>) (rjo being the 
suspending fluid viscosity and <f> the- volume fraction). 
These theories predict a negative value for the coefficient 
K for pushers at low shear rates, meaning the suspension 
can exhibit a lower viscosity than the suspending fluid. 
The theoretical assessment of shear viscosity relies on an 
assumed statistical representation of the orientations of 
the bacteria, captured by a Fokkcr-Plank equation and 
a kinematic model for the swimming trajectories [25I |26| . 

Despite the large number of theoretical studies, few 
experiments have been conducted. With bacillus sub- 
tilis (pushers) trapped in a soap film Sokolov et al.(l2l] 
have shown that a vorticity decay rate could be asso- 
ciated with a strong decrease of shear viscosity in the 
rcsence of bacteria. For algca (pullers), Rafai et al. 
have shown that the effective viscosity measured in 
a classical rhcometer is larger than the viscosity of the 
corresponding dead (passive) suspension. However, no 
measurements of the viscosity of a dilute suspension of 
pushers under controlled shear conditions exist to date. 
This is mainly due to the fact that one has to assess 
small viscosities near the viscosity of water at very small 
shear rates to probe the theoretical predictions. These 
parameters are typically outside of the resolution of a 
classical rotational rhcometer and have thus made these 
measurements inaccessible. In this letter, we present the 
first measurements of the shear viscosity of a suspension 
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of pushers using a microfluidic device to overcome these 
difficulties and obtain the relative viscosity of an active 
suspension for a large range of shear rates and bacteria 
concentrations. Our set-up also allows direct visualiza- 
tion of the flow as well as the spatial distributions of 
bacteria in the flow. 

The wild type E. coli W used here are prepared follow- 
ing the experimental procedure described in Ref . 0, [H| . 
The strain is grown overnight in rich medium (LB). After 
washing, it is transferred into MMA, a motility medium 
supplemented with K-acetate (0.34 mM) and polyvinyl 
pyrolidone (PVP: 0.005%). The sample is then incubated 
for at least one hour. To avoid bacteria sedimentation, 
Percoll is mixed with MMAP 1 vol/ 1 vol (isodense con- 
ditions). We verified that under these conditions, the 
suspending fluid is Newtonian (viscosity rj = 1.28 x 10~ 3 
Pas at 22° C). All experiments are performed at a fixed 
temperature T = 25°C. 

To obtain the shear viscosity, we adapted a microfluidic 
device [l?} to compare the Newtonian viscosities of two 
liquids. The device is a Y-shape Hele-Shaw cell of height 
h, such that each branch receives a different fluid (resp 
fluid of viscosity 770 and fluid 1 of viscosity 771 (see figH]) . 
Both flows are driven at an identical flow rate Q. In a 
Hele-Shaw approximation, the velocity profile in the z di- 
rection can be described by a parabolic flow profile and 
the dominant shear rate occurs in the z direction. Under 
these conditions and for a viscosity ratio near one, the 
position of the interface between both fluids at steady 
state i.e. di/do is determined by the viscosity ratio 771/770 



271 ]. Guillot et al. [27( have also used this approach for 



non-Newtonian fluids, where it corresponds to the mea- 
surement of an apparent viscosity. Here we will follow the 
same approach and will verify our method subscquanty. 
The suspension of bacteria is flown into one arm and the 
suspending Newtonian fluid into the other arm. The in- 
terface position is then measured at various flow rates Q. 
The experimental data is presented as a function of the 
maximum shear rate obtained by assuming a parabolic 
flow profile 7 M = (6Q)/(h 2 di) where d\ is the lateral 
width occupied by the suspension. The relative viscosity 
T] r = is then : 
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This micro-fluidic device has the advantage of measur- 
ing a viscosity ratio and provides very good resolution 
of the suspension viscosity independent of its absolute 
value or the applied shear rate. Special care has to be 
taken that both branches are fed at the same constant 
flow rate Q. We use a very high precision two-syringe 
pump from nemeSYS and a precision syringe (Hamilton 
Gastight 1805RN) of a very small volume (50/iL) allow- 
ing to impose identical and very small flow rates (down 
to Q = 0.5 nL/scc) on both arms. Finally, the posi- 
tion of the sample region (indicated by the red rectangle 



on figure []} is chosen in such a way as to be in steady 
state conditions (reached at a distance of approximately 
600 fj,m from the junction of the inlet channels) while 
avoiding significant mixing. 

The Y-shape channel was fabricated completely from 
PDMS using a soft-lithography technique. The channel 
thickness is h = 100 fim. The main channel width is 
w = 600 fim and the two inlet branches widths are w/2. 
Inlets are connected by 500 /im diameter tubes to the 
two-syringe pump. The total length of the main chan- 
nel is 40mm. Suspensions were prepared (see Q) with 
a number of bacteria per unit volume n in the range 
1.9 lO 12 ^- 1 < n < 25.6 10 12 £ _1 . The volume fraction 
is estimated using the space occupied by the body of 
each bacteria Vb = 1 fim 3 such that </> = n Vb, yielding a 
range of 0.19 < <f> < 2.56%. Note that other definitions 
of the volume fraction, for example based on the length 
of the body or including the length of the flagella would 
lead to significantly higher volume fractions. The flowing 
suspension was observed using an inverted microscope 
(Zeiss-Observer, Zl) connected to a digital camera Pix- 
eLINK PL-A741-E (1280 X 102Apix 2 ) capturing videos 
at a frame rate of 22 images/s. Low magnification (2.5x) 
allowed an extended view of the channel (sec FigQ]). 
During an experiment we begin with an initial flow rate 
of Q = 0.5 1 nL/sec, and then the input flow is increased 
step by step up to Q = 100 nL/sec. Note that we have 
verified for all experiments that identical results are ob- 
tained when subsequently decreasing the flow rate. In fig 
[TJ the shape of the interface obtained by averaging over 
120 successive images is displayed for the measurement 
area. To quantitatively determine the interface position, 
we measure the mean bacteria density across the chan- 
nel width n(z) using the logarithm of the profile of the 
average intensity C(y) = \n((I)(y)/ {I)o as a measure of 
the concentration, with (I)q is the mean intensity in the 
absence of bacteria. The average is calculated in the x- 
direction over a distance of 600 fim (see red rectangle 
in figH]). The concentration C(y) is fitted with an error 
function erf(?/) to obtain the interface position yi. Once 
the interface position is determined, we extract the rela- 
tive viscosity according to equation |T]) and associate it 
with j M . 

Experimental observations are presented in fig. [2]Ja) dis- 
playing averaged images for flow rates of Q = 0.5 nL/sec, 
Q = 2 nL/scc and Q — 10 nL/sec for a suspension 
at a volume fraction of <f> = 0.35%. When increas- 
ing the flow rate the interface position changes from a 
value above the mid-position to a value below the mid- 
position indicating a change in the suspension viscosity 
from greater then the viscosity of the suspending fluid to 
smaller. Quantitative measurements are given in figHJe) 
showing the relative viscosity r\ r of the suspension as 
a function of the maximum shear rate for various con- 
centrations 6 = 0.19, 0.75, 2.3 and 2.48 %. The er- 
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FIG. 2: Rheology curves. Fig (a) are 3 averaged pictures of 
the bilaminar flow for flow rates of Q — 0.5 nL/s, Q — 2 nL/s 
and Q = 10 nL/s respectively for a (f> — 0.35% bacteria sus- 
pension. These images are averaged over 120 images at 22im/s 
and observed with a 2.5x magnification using phase contrast. 
The dotted line represents the mid-position in the channel 
width. The scale bar corresponds to 200 um. (b) Relative 
viscosity versus mean shear rate at several volume fractions 
(c) Relative viscosity versus volume fraction at several shear 
rates. The errors bars are estimated using the detection error 
of the interface position. 



ror bars are estimated using the uncertainty in deter- 
mining the interface position and in marked contrast to 
classical rheology, the error does not increase for lower 
shear rates, once again demonstrating the strength of our 
micro-rheometer. For the lowest concentration the vis- 
cosity of the Newtonian suspending fluid is recovered for 
all shear rates validating our rheological device. For all 
other concentrations, the curves display the same qual- 
itative behavior. At low shear rates, we observe shear 
thickening and at higher shear rates, shear thinning oc- 
curs. For all of these cases, the shear thinning and shear 
thickening character of the active suspensions are weak 



and power law indices n (using 77 = Kj^ n 1 ^) close to one 
are found for both regimes for all concentrations. At all 
concentrations tested, a maximum viscosity is observed 
at a value around j M = 20s~ 1 . The non-monotonic be- 
havior as a function of shear-rate is in agreement with 
the results of Saintillan [24} obtained for slender bacte- 
ria. Note that the low shear rate plateau predicted by the 
theory could not be observed in our experiments due to 
practical limitations of our device (minimum flow rates 
0.5 nL/sec). At low volume fractions and low shear rates, 
we observe a decrease of rj r below one which agrees with 
the theoretical predictions of Haines et al. [22| and Sain- 
tillan [24| of this striking effect. Figure [He) shows the 
relative viscosity as a function of the volume fraction <j> at 
various shear rates 7m = 1, 2, 4 and 20 s _1 . A decrease 
in viscosity is observed for the small shear rates, followed 
by a sharp increase of viscosity for all shear rates, corre- 
sponding to a semi-dilute regime. In our case this regime 
is observed for concentrations above approximately 1%. 
Similar behavior was also observed by Sokolov et al.[l2[ 
using vortex decay in a suspension of Bacilus subtilis con- 
fined in a soap film and has been predicted by Ryan et 
al. [23| in their simulations. 

So far we have measured the shear viscosity of the suspen- 
sions as a function of the maximum shear rate assuming 
a parabolic velocity profile. One of the advantages of our 
microfluidic device is that we can access the local veloc- 
ity and bacteria concentration profiles as a function of 
the channel height directly. Bacteria moving in the flow 
were visualized with a high magnification objective (40x, 
phase contrast) allowing the position of the bacteria to 
be monitored at various heights z (field depth 3 fim). 
Videos were taken using a highspeed camera (Photron 
FastCam SA3, resolution 1024xl024pixels, shutter speed 
(l/500s), frame rate l/50s for Q = 0.5, lnX/s and 
l/500s for Q = lOnL/s). To monitor the flow ve- 
locity, we suspended a very low concentration of 2/im 
density-matched latex beads as passive tracers (density 
1.03g/cm 3 ). Here we present selected results for the 
volume fraction <f> = 0.75 % at which the reduction of 
the relative viscosity below 1 is observed at three flow 
rates {Q = 0.5,1,10 nL/sec corresponding respectively 
to maximum shear rates jm = 1, 2, 20 s _1 . Videos were 
taken in a region situated in the middle between the in- 
terface and the side wall (see blue frame of fig. [T|). To 
reduce the detection noise, a bacterium is retained only if 
it is detected on at least two consecutive frames. The flow 
velocities were computed for particles moving in the focal 
plane. In Fig. [3^,) we see that the concentration profiles 
are similar in shape to those published previously by 15 1 
in the absence of flow, i.e. a quasi-constant density and 
a strong density increase within lOfim of the wall, hence 
probing a trapping effect that persists over the range of 
shear rates explored. The calculated mean concentration 
is represented by the dashed line in fig. [5Ji and is in agree- 
ment with the measured bulk concentration. It is impor- 
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tant to note that the concentration profiles are identical 
for the various shear rates. The change in viscosity ob- 
served is therefore not due to a change in concentration 
within the microchannel. The velocity profiles V x (z) are 
displayed in Fig.([5]3). Within the precision of our exper- 
imental setup, a deviation from a parabolic velocity pro- 
file is not observed justifying the approximation used for 
our analysis. Note also that this finding is in agreement 
with the weak shear thinning or shear thickening char- 
acter of our suspensions, that predict only small changes 
in the velocity profile, not detectable in our experiments. 
For each flow rate the curves presented in fig. |3p) agree 
within 3% with the theoretical predictions for a rectan- 
gular cell of the aspect ratio of (1/6) fij^ . 
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FIG. 3: Bulk profile measurements for a volume fraction cj> = 
0.75 % and shear rates 7m = 1, 2, 20 s . (a) Bacteria density 
profile n(z). The horizontal dashed line represents the mean 

density: n = 7.5bact/L. (b) Velocity profile of passive tracers 
VJz). 



In conclusion, using a micro-fluidic device, we have 
measured for the first time the effective shear viscosity 
of a suspension of "pushers" ( E-coli bacteria) over a 
large range of shear rates(l — 200 s _1 ) under controlled 
shear conditions in the dilute and semi-dilute regime. 
We confirmed an important prediction for the rheology 
of pushers: the active viscosity can be lower than the 
viscosity of the suspending fluid at low shear rates 
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24j. In the dilute as well as in the semi-dilute 
regime, we observed a shear-thickening behavior at 
lower shear rates followed by a shear-thinning regime at 
higher shear rates. The viscosity maximum is observed 



at a shear rate on the order of the inverse of the time a 
bacterium needs to swim over a distance of its own size 
and this value seems independent of the bacteria concen- 
tration. These results are consistent with the theoretical 
calculation of Saintillan 
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Direct measurements of the 
bacteria concentrations and flow velocities confirm that 
the bacteria distribution remains homogeneous and the 
parabolic flow profile is not significantly modified by the 
weak non-Newtonian character of the active fluid. This 
last result validates the use of the Newtonian approxima- 
tions for the extraction of our viscosity measurements. 
In the semi-dilute regime (here for volume fractions 
greater than 1%), we observe a strong increase of the 
viscosity consistent with numerical simulations by Ryan 
et al.[23j. Our results represent the first experimental 
validation of the non-Newtonian rheology of an active 
suspension of pushers under controlled shear conditions. 
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